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We discuss some recent results on the statistics of the Coulomb Blockade in disordered quan- 
tum dots containing spinless interacting fermions using the self-consistent Hartree-Fock ap- 
proximation. We concentrate on the regime r s > 1, with finite dimensionless conductance g. 
We present significantly different results for the cases of a Coulomb and a nearest-neighbour 
bare interaction. We show that neglecting rearrangements when the particle number is 
changed (Koopmans' approximation) can lead to large errors, and show how the spectral 
statistics of the self-consistent single particle orbitals evolve with the interaction strength. 



1 Introduction 

The Coulomb Blockade to tunnelling in quantum dots, where the charge in the dot is a good 
quantum number, provides an experimentally accessible measure of the difference between 
ground state energies of successive total particle number mm. More specifically, the ac- 
cessible quantity is the second difference of ground state energies Eq[N, V g ): 

A 2 (N) = E G (N + 1, V G ) - E G (N, V G ) - E G (N, V G ) + E G {N - 1, V G ) , (1) 

where V G parametrises the gate voltage environment at each resonant tunnelling point. Within 
the constant interaction (CI) model the ground state energy is simply the sum of filled single 
particle energies £ (n) plus N(N — l)e 2 /2C, where e 2 /C is the constant interaction, and C oc L _1 
is the geometric capacitance. Taking V G = V G , the peak spacing trivially reduces to 

A 2 (N) = £{N + l)-£{N) + e 2 /C (2) 

and so, in the diffusive regime, displays Wigner-Dyson (WD) statistics! shifted by e 2 /C, up to 
corrections in one over the dimensionless conductance, g. 

Recent experiments on quantum dotsBBiH have shown that whilst the mean peak spacings 
are apparently well described by the CI model, the fluctuations are not described by Wigner- 
Dyson statistics. It is found that the distribution of A 2 is roughly Gaussian EJtra, with broader 
non-Gaussian tails seen in Ref. 3. In Refs. 1,2,4 the variance of the fluctuations was found to 
be considerably larger than that given by the WD distribution. Exact numerical calculations on 
very small systems □ also exhibit Gaussian fluctuations on a scale much larger than A oc L~ 2 , 
and it was suggested that the typical size of the fluctuations scale universally with the average 
spacing (A 2 ) 2 e 2 /C. 

In an attempt to bridge the gap between theory a and experiment, Blanter et aM evaluated 
the fluctuations within a Hartree-Fock framework, neglecting effects due to the change in gate 
voltage V G to V G . To this end they applied the Random Phase Approximation (RPA) to gener- 
ate the screened interaction in the confined geometry, and assumed that all HF level spacings, 
except for the Coulomb gap, are described by WD statistics. Implicitly assuming Koopmans' 
theorem t3 to be valid (which neglects single-particle rearrangements when the particle number 
distribution is altered), and using wavefunction statistics established for non- interacting elec- 
trons in a random potential, they calculated the fluctuations of A 2 beyond the CI model. These 
additional fluctuations were found to be parametrically small (in 1/g) and proportional to A. 
Hence, the total fluctuations in A 2 were found to be proportional to A. The disagreement 



between experiment and theory therefore remained. The analysis of Ref. 9 is consistent in the 
limits g > 1, r s C 1. The parameter r s characterises the relative importance of interactions 
in the electronic system, and is defined as the mean electron separation in units of the effective 
Bohr radius. Typical values of r s ~ 0(1) have been reported in experiment §@. 

To understand why theory and exact numerics are in disagreement, we consider an effective 
single particle model as in Ref. 9, but generate the single particle orbitals self-consistently 
using the Hartree-Fock approximation (SCHF). This model allows for the rearrangements when 
particles are added or removed which was neglected in Ref. 9, and does not require any conjecture 
about the wavefunction correlations (the two issues are of course intimately linked). We are 
also able to consider much larger samples than is feasible by exact diagonalisation 0. This 
approximation has also been seen to be quite good for the calculation of persistent currents in 
similar systemslij. We show that fluctuations large compared to the single particle level spacing 
can arise without recourse to varying the sample shape, size or gate-to-dot coupling, supposing 
these to be additional effects. We further demonstrate that approximating the addition spectrum 
spacings by applying Koopmans' theorem can lead to large errors in the calculation of the spacing 
statistics. A more detailed presentation of these results can be found elsewhere liJtJ. 



2 The Model 



We consider the following tight binding Hamiltonian for spinless fermions 

H = Y, Wicfd - t 4+ri C i + 2" 12 M ijij C t C j C i Ci ( 3 ) 
i i,7] ij 

where i is the site index, rj describes the set of nearest neighbours, Wi is the random on site 
energy in the range [-W/2, W/2], and t the hopping matrix element, henceforth taken as unity. 
We study separately, both a Coulomb interaction potential, 

Mijij = l/|rj ~ r j\ , (4) 

and a short range potential, 

M-ijij = <5« , « —i— 77 , (5) 

there is no need to add a constant interaction as it does not affect the physics. 

We consider a 2D system with periodic boundary conditions, and choose to define 

|rj - rj\ 2 = [L x sm 2 (Trn x /L x ) + L 2 $m 2 (irn y / 'L y )]/ir 2 , (6) 

where (n x ,n y ) = rj — rj, and L x , L y are the sample dimensions. 

All distances are measured in units of the lattice spacing a; the physical parameters are 
therefore U = e 2 /a, t = h 2 /2ma 2 . The standard definition for r s is given, for low filling, by 
r s = V I (t\/4:Trv), where v = N/A is the filling factor on the tight binding lattice with A sites. 
The dimensionless conductance g can be approximated, again for low filling, using the Born 
approximation. We find g = Q&irv{t/W) 2 , which is valid for A, N ^> g ^> 1. 



3 The SCHF Approximation 

There are a few aspects of this approximation scheme which warrant discussion before presenting 
the results. Firstly we should stress that this scheme is non-perturbative, but uncontrolled. To 
ensure that the converged solutions obtained were at least local minima, and hence the RPA 
excitations stableB, we diagonalised the total energy curvature matrix0 in some representative 
samples. The calculational effort involved is too heavy to do this for all cases. 



cuiation of the RPA screened 
fl It was found in Ref. 9 that 



The issue of which bare potential (f|,||) should be used is subtle. If screening is generated 
externally (by close metallic plates for example), it is reasonable to insert a short-ranged inter- 
action for all interaction strengths. It seemsEJ that this model may then be useful to model the 
quantum dots of Refs. 5,6. It is not clear however that this form should be taken to represent 
screening by high-energy particle hole pairs (the RPA approximation). The reason being that 
the self-consistency process automatically incorporates much of the physics included in the RPA 
process. In other words, it generates many of the diagrams included in the RPA calculation, 
and so leads to double counting. This problem does not arise if one applies the Coulomb bare 
potential. However, we know from basic considerations that the exchange (Fock) interaction 
should be screened, and the SCHF process is unable to achieve this. In clean systems this can 
lead to large errors, here however, we have verified that the typical fluctuations of the exchange 
contribution are smaller than those of the Hartree contribution over the parameter range con- 
sidered. This suggests that the error made in not screening the exchange term correctly is not 
overly important. We expect then, that the results for a long ranged interaction should be 
representative of the exact numerical calculations. 

The simplest way to interpret the RPA calculation, is the (average) response of the electron 
gas to the insertion of a test charge e. Since the excess charge cannot be pushed to infinity in 
a bounded dot, it is perhaps not surprising that an explicit a 
interaction shows that a charge e accumulates on the boundary 
such boundary effects could dominate the interaction induced fluctuations of A2. In any event, 
the effect of the boundaries is geometry specific; we remove it by considering a torus geometry. 

Since the SCHF approximation attempts to generate static screening by rearranging the 
effective-single-particle orbitals, there must be a different interpretation between the SCHF 
approximation with a Coulomb bare potential, and HF with RPA. In the latter, most charge 
rearrangements are made via the RPA summation, and as such, are separated from the effective- 
single particles. In our time-independent self-consistent theory, it is not possible to make such a 
separation, and indeed the job of charge rearrangement is carried purely by the effective-single- 
particle orbitals. Clearly then, the meaning of the effective-single-particle is not the same in the 
two cases. 

Since the validity of Koopmans' theorem is analysed here, we explain our motivation for 
investigating an effect that is second order effect in r s . At sufficiently small r s , where non- 
interacting wavefunction correlations can be employed, the effects of rearrangements are neg- 
ligible. It is often assumed that they are small for the further reason that adding a particle 
generates only O(N) corrections to the existing single-particle wavefunctions. We stress how- 
ever, that when evaluating the ground-state energy, one must sum over N such corrections (as 
well as iV 2 of 0(N 2 )), and so the cumulative effect can be large. When r s is sufficiently large as 
to affect the wavefunction correlations, it is further plausible, that the 0(r 2 ) corrections to the 
wavefunctions correlate with the unperturbed wavefunctions to generate typical fluctuations of 
order r s . 



4 Nearest-Neighbour Interactions 



In Ref. 11 we show that the distribution of level spacings for both occupied and unoccupied 
states evolve with r s from a Wigner-Dyson distribution towards a Poisson-like distribution. 
This is more pronounced for the unoccupied states, which only have a meaning if Koopmans' 
approximation is invoked. This tendency is rather weak however, and for r s ~ 0(1), the 
distribution is almost indistinguishable from Wigner-Dyson. On the other hand, we found 
that the self-consistently calculated gap (A2) evolved comparatively rapidly towards a Gaussian 



a We note that if the dot is supposed to be charge neutral, it may be more consistent if the electron gas were 
to see an excess background charge of — e, plus the test electron in question 



distribution, and was essentially Gaussian for r s > 2. The main difficulty with observing the 
typical interaction induced fluctuations for r s ;$ 1 was that they were small compared to A. In 
this regime we confirmed that they scaled like r s A for fixed g, in agreement with Ref. 9. Indeed, 
if the wavefunctions remain ergodic, it is clear that this scaling must hold. It was too difficult 
to identify the disorder scaling with any certainty. In this regime, applying Koopmans' theorem 
did of course generate errors, but they were found to be relatively small: the mean error scaled 
like rl A, and the error in the typical fluctuations was too small to analyse with confidence. 

When r s > 1, we found a dramatic development in the SCHF ground state, and a correspond- 
ingly dramatic growth in the error incurred through Koopmans' approximation. We also noted 
that the system size scaling of the mean Koopmans' error only settled to the r^A dependence 
for the larger systems considered (^4 > 50 lattice sites); well beyond the current limitations of 
exact calculations. Whilst the single-particle wavefunctions remained essentially ergodic, incipi- 
ent charge density modulations (CDMs) led to an unusual, and highly non-universal fluctuation 
behaviour. In this regime, deviations in (A2) from the CI model prediction were observed, and 
identified to be 0(U/A). They were shown to be a direct consequence of the incipient CDMs. 
These CDMs are analysed in much greater detail in Ref. 13. 

In conclusion, we find that for weak interactions, the typical fluctuations 5Ao ~ f s A for 



incipient charge density modulations cause a significant change in the statistics of A2. It is 
suggested in Ref. 13, that such effects may be relevant to the unusual addition spectra observed 
in the capacitance experiments of Refs. 5,6. 

5 Coulomb Interactions 

As for the nearest-neighbour interaction case, the distribution of level spacings for both occupied 
and unoccupied states were investigated in Ref. 11. They were both found to remain very close 
to the WD distribution for all interaction strengths considered (up to r s ~ 5). The gap (A2) 
distribution was found to evolve towards a Gaussian distribution much more rapidly than changes 
in the above level spacings, but the scale on which this happened depended on the size of the 
sample. 

The mean error generated by Koopmans' approximation was found to tend towards a simple 
quadratic dependence on r s as the system size was increased to the maximum size considered 
of 144 sites and 37 particles. As for the nearest neighbour case, the system size scaling changed 
around A ~ 50 lattice sites. The mean Koopmans' error therefore diverges with system size when 
compared to both (A2) and A. The Koopmans' approximation error in the typical fluctuations 
was found to become very large (0(<5A2)), no attempt was made to generate an empirical formula 
to describe it. 

The divergent relative error of Koopmans' approximation to the mean charging energy was 
one of the more striking results of this work. Whilst the CI model worked remarkably well, a 
supposed improvement on it (Koopmans' theorem about a self-consistent basis) was found to 
be worse, and moreover, could only be made meaningful for r s < 0{1/L). This contradicts 
the expectation that Koopmans' theorem should become a good approximation for r s <C 1, 
regardless of the system size. In the CI model this error vanishes identically, while in Blanter et 
alu, the terms which generate it were implicitly assumed to vanish. 

Turning to the main object of investigation, the SCHF approximation was found to yield 
typical fluctuations in the charging energy of 5A2 ~ 0.52A + a(A2) / \[g + 0{r 2 s ), where a is 
a constant. A similar result was also found independently by Bonci and Berkovits 0. Whilst 
this result disagrees with the claim in Ref. 1 that the fluctuations are universal, the system 
size scaling is consistent. The 0(rf) contribution arises due to interaction induced correlations 
between the single-particle wavefunctions, and eventually dominates the fluctuations as r s is 



fixed g > 1, in broad agreement with Blanter 




increased. However, the validity of the SCHF approximation becomes more questionable as r s 
is increased, and it is not clear to what extent the results are valid in this strong interaction 
regime. One might expect that the combined fluctuations of the single particle orbitals and 
the RPA-bosons in the HF with RPA model are similar to the fluctuations generated in the 
Coulombic SCHF model. Since we do not find agreement with Blanter et alu, this appears not 
to be the case. A direct calculation of fluctuations in the polar isability 01 does not, account 
for this discrepancy. It would appear that either screening is not generated in our system (or 
indeed in the exact calculations on very small samples), or that it is not correct to separate out 
the RPA-bosons from the fluctuation effects in the the usual way. 

In summary, we found that Koopmans theorem fails dramatically in the case of a Coulomb 
bare interaction. We find that the typical interaction induced fluctuations are proportional to 
(A2) to leading order in r s , but they are not universal, rather they vanish with the dimensionless 
conductance like 1/y/g- Other related work has been independently carried out by Levit and 
OrgadHa, Cohen, Richter and Berkovits Ej, and Bonci and BerkovitsliS. 
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